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Abstract 

The reaction pp — > pf( y 7T + 7i^'n- + 7i^)p s has been studied at 450 GeV/c in an experiment 
designed to search for gluonic states. A spin analysis has been performed and the cLPt filter 
applied. In addition to the well known /i(1285) there is evidence for two J PC = 2~ + states 
called the r] 2 (1620) and rj 2 (l875) and a broad scalar called the /o(2000). The production of 
these states as a function of the dPx kinematical filter shows the behaviour expected for qq 
states. In contrast, there is evidence for two states at 1.45 GeV and at 1.9 GeV which do not 
show the behaviour observed for qq states. 



Dedicated to the memory of our colleague and friend Yuri Prokoshkin 



Submitted to Physics Letters 



* Deceased. 

1 LAPP-IN2P3, Annecy, France. 

2 Athens University, Nuclear Physics Department, Athens, Greece. 

3 Bergen University, Bergen, Norway. 

4 School of Physics and Astronomy, University of Birmingham, Birmingham, U.K. 

5 CERN - European Organization for Nuclear Research, Geneva, Switzerland. 



IHEP, Protvino, Russia. 
IISN, Belgium. 
JINR, Dubna, Russia. 

High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305, Japan. 
Faculty of Engineering, Miyazaki University, Miyazaki, Japan. 
RCNP, Osaka University, Osaka, Japan. 
Oslo University, Oslo, Norway. 

Faculty of Science, Tohoku University, Aoba-ku, Sendai 980, Japan. 
Faculty of Science, Yamagata University, Yamagata 990, Japan. 



The WA76 and WA91 collaborations have studied the centrally produced 7r + 7r~7r + 7r~ final 
state in the reaction 

PP — > Pf(7T + 7T~7T + 7T~)p s (1) 

at 85 0, 300 0, and 450 GeV/c ||. The subscripts f and s indicate the fastest and slowest 
particles in the laboratory respectively In addition to the /i (1285), which was observed at 
all energies, the WA76 collaboration observed two new states at 1.45 and 1.9 GeV in their 
300 GeV/c data 0. In contrast, no clear evidence was seen for these states in the 85 GeV/c 
data of the same experiment []J . The increase in cross section with increased incident energy |J 
is consistent with the formation of these states via a double Pomeron exchange mechanism, 
which is predicted to be a source of gluonic states |§]. 

The WA91 experiment ||, which studied reaction (1) at 450 GeV/c, confirmed the existence 
of these states and, from a spin analysis, showed that the peak at 1.45 GeV had I G J PC = + ++ 
and that the peak at 1.9 GeV had I G J PC = + 2 ++ ; hence they were called the /o(1450) and 
/2(1900) respectively. In a reanalysis of this data ||, the WA91 collaboration showed that it 
was possible to describe the /o(1450) as being due to the interference between the /o(1300) and 
the / (1500). 

This paper presents new results from the WA102 experiment which is a continuation of the 
WA76 and WA91 programme and has more than a factor of five increase in statistics. This 
enables a detailed spin analysis to be performed and also a study as a function of dPp, which is 
the difference in the transverse momentum vectors of the two exchanged particles and has 
been proposed as a glueball-gg filter @ . 

The data come from experiment WA102 which has been performed using the CERN Omega 
Spectrometer. The layout of the Omega Spectrometer used in this run is similar to that 
described in ref . || with the replacement of the OLG A calorimeter by GAMS 4000 P . 

Reaction ([]]) has been isolated from the sample of events having six outgoing tracks by first 
imposing the following cuts on the components of the missing momentum: |missing P x \ < 14.0 
GeV/c, |missing P y \ < 0.12 GeV/c and |missing P z \ < 0.08 GeV/c, where the x axis is along the 
beam direction. A correlation between pulse-height and momentum obtained from a system of 
scintillation counters was used to ensure that the slow particle was a proton. 

The quantity A, defined as A = MM 2 (pfp s ) — M 2 (7t + 7t~7t + tt~), was then calculated for 
each event and a cut of |A| < 3.0 (GeV) 2 was used to select the 7r + ix~ix+ix~ channel. Events 
containing a fast A ++ (1232) were removed if M(pfii + ) < 1.3 GeV, which left 1 167 089 centrally 
produced events. 

Fig. |TJa shows the acceptance corrected 7r + 7r~7r + 7r~ effective mass spectrum renormalised to 
the total number of observed events. The mass spectrum is very similar to that observed by 
experiments WA76 and WA91, namely, a clear peak at 1.28 GeV associated with the /1 (1285), 
a peak at 1.45 GeV called the /o(1450), which is possibly due to an interference effect between 
the /o(1300) and /o(1500), and a broad enhancement at 1.9 GeV called the /2(1900). The 
mass spectrum has been fitted using three Breit-Wigners representing the /1 (1285), the 1.45 
GeV peak and the ^(1900) plus a background of the form a(m — mth) b exp(—cm — dm 2 ), 
where m is the 7r + 7r~7r + 7r _ mass, m t h is the 7r + 7r~7r + 7r~ threshold mass and a, b, c, d are fit 
parameters. Reflections from the r]7i + 7i^ decay of the rj and /i(1285) give small enhancements 
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in the 7r + 7r~7r + 7r~ mass spectrum in the 0.8 and 1.1 GeV regions due to a slow ir° from the 
decay of an r\ falling within the missing momentum cuts. In order to get a correct description of 
the 7T + 7r~7r + 7r~ mass spectrum two histograms representing a Monte Carlo simulation of each 
reflection have been included in the fit. The Breit-Wigner describing the /i(1285) has been 
convoluted with a Gaussian with a = 12 MeV representing the experimental resolution in this 
mass region. The masses and widths determined from the fit are 

M x = 1281 ± 1 MeV I\ = 24 ± 3 MeV 
M 2 = 1445 ± 4 MeV T 2 = 95 ± 15 MeV 
M 3 = 1920 ± 20 MeV T 3 = 450 ± 60 MeV. 

Close and Kirk |7[ have proposed that when the centrally produced system is analysed as 
a function of the parameter dPx, which is the difference in the transverse momentum vectors 
of the two exchange particles [EJ, states with large (small) internal angular momentum will be 
enhanced at large (small) dPx- A study of the 7r + 7r~7r + 7r~ mass spectrum as a function of dPr 
is presented in fig. |b), c) and d) for dP T < 0.2 GeV, 0.2 < dP T < 0.5 GeV and dP T > 0.5 GeV 
respectively. A dramatic effect is observed; the /i(1285) signal has virtually disappeared at low 
dP T whereas the / (1450) and / 2 (1900) signals remain. 

A fit to figs. [TJb) , c) and d) has been performed (not shown) using the resonance parameters 
fixed from the fit to the total mass spectrum. A poor quality fit results which indicates that 
the underlying resonance structure may be more complicated. In order to try to understand 
the resonances that are present in the system an extraction of the partial waves is required. 

A spin-parity analysis of the 7r + 7r~7r + 7r~ channel has been performed using an isobar model ||. 
Assuming that only angular momenta up to 2 contribute, the intermediate states considered 
are 



aa, a(7l + 7T )swave, a(7T + 7l ) 

Pwavei 0"(7T^7T ^Dwavei O P i 

P°P°, P°(^ + ^ )swave, P°(7r + 7r- ) Pwave , p° (n + 71^ ) Dwave , 

ai (1260)7T, a 2 (1320)7r, / 2 (1270)<r, / 2 (1270)(7r+7r-) 5 ^ e , 

/ 2 (1270)p°, / 2 (1270)/ 2 (1270) 



where a stands for the low mass tttt S-wave amplitude squared. Two parameterisations of the 
a have been tried, that of Au, Morgan and Pennington JTIJ and that of Zou and Bugg [II 



The amplitudes have been calculated in the spin-orbit (LS) scheme using spherical harmonics. 

In order to perform a spin parity analysis the log likelihood function, Cj = l°g -P? W 5 i s 
defined by combining the probabilities of all events in 40 MeV 7r + 7r~7r + 7r~ mass bins from 1.02 
to 2.82 GeV. The incoherent sum of various event fractions cij is calculated so as to include 
more than one wave in the fit, 

^ = E lo g[E^) + (i-E%)) ( 2 ) 

i \ J i I 

where the term (1 — J2j %) represents the phase space background. The negative log likelihood 
function (— C) is then minimised using MINUIT Coherence between different J p states 



and between different isobar amplitudes of a given J have been neglected in the fit. Different 
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combinations of waves and isobars have been tried and insignificant contributions have been 
removed from the final fit. As was already mentioned in the WA91 publication ||, in the 
analysis of the /i (1285) and the / (1450) peak there is little change in the result if the pp 
amplitude is used instead of the p(inr)p wave . Using Monte Carlo simulations it has been found 
that the feed-through from one spin parity to another is negligible and that the peaks in the 
spin analysis can not be produced by phase space or acceptance effects. 

The results of the best fit are shown in fig. |2| for the total data sample and in fig. |3] for 
three ranges of dPr- The /i(1285) is clearly seen in the J p = 1 + pp. Superimposed on the 
J p = 1 + pp wave for the total sample is the Breit Wigner convoluted with a Gaussian used 
to describe the /i(1285) in the fit to the mass spectrum. As can be seen the /i(1285) is well 
represented in both size and shape. From fig. |3] it can be seen that the /i(1285) dies away for 
dP T < 0.2 GeV. 

The J p = + pp distribution in fig. ||] shows a peak at 1.45 GeV together with a broad 
enhancement around 2 GeV. In order to test if the peak at 1.45 GeV is due to ao rather than 
pp the change in log likelihood in the four 40 MeV bins around 1.45 GeV has been calculated 
by replacing the J p = + pp amplitude by the J p = + aa amplitude. The likelihood decreases 
by AC = 9000 corresponding to n — \/2A£ = 134 standard deviations. This raises a possible 
problem in interpreting the peak at 1.45 GeV with the /o(1500) since the analysis of the 4n 
channel in pp by Crystal Barrel |L3| and in radiative J /if) decay [0 shows that the /o(1500) 



decays via aa. A possible explanation could be that the interference between the /o(1300) and 
/o(1500) produces the observed pp decay. This question should be answered in the future by a 
study of centrally produced 47r° events. 

From fig. [3] it can be seen that the peak in the J p = + pp wave around 1.45 GeV remains 
for dPx < 0.2 GeV while the J p = + enhancement at 2.0 GeV becomes less important: which 
shows that the dPr effect is not simply a J p filter. A fit has been performed to the J p = + pp 
amplitude in fig. Q in the 1.45 GeV region using a single Breit- Wigner and gives a mass of 
1445 ± 4 MeV and a width of 200 ± 20 MeV. The mass is compatible with the fit to the total 
mass spectrum although the width is broader; it is nonetheless compatible with a fit performed 
to the mass spectrum for dPx < 0.2 GeV. A similar fit has been performed to the J p = + pp 
distributions shown in fig. |3|. The mass and width found in each interval are constant within the 
errors of the fit. This would indicate that either the peak is due to a single resonance or that 
if it is due to interference between the /o(1300) and /o(1500) then both states have a similar 
dPx dependence. 

In order to test the interference hypothesis a fit has first been performed to the total 
J p = + pp distribution in fig. |2] using a K matrix formalism [|15| including poles to describe 



the interference between the /o(1300), the /o(1500) and a possible state at 2 GeV. The result 
of the fit is superimposed on the J p = + pp distribution shown in fig. [2] and well describes the 
data. The resulting resonance parameters are given in table [|. The parameters of the /o(1300) 



and /o(1500) are very similar to those found by Crystal Barrel ]I6 



As can be seen from fig. § neither the J p = 2 + a 2 (1320)7r nor the J p = 2 + J * 2 (1270) (nil) swave 
alone can describe the /a(1900) peak observed in the mass spectrum. However, the sum of 
the two waves accounts for most of the /2(1900) signal. The a2(1320)7r distribution and the 
/2 (12 70 )(vr7r) swave distribution peak at different masses which suggests that the / 2 (1900) could 
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be composed of two J PC = 2 ++ resonances but neither distribution is well represented by a 
Breit-Wigner. However, the sum of the J p = 2+ a 2 (1320)vr and J p = 2+ f 2 (l270)(7nr) Swave 
waves has been fitted with a single Breit-Wigner, giving a better description, with M = 1960 
± 30 MeV and T = 460 ± 40 MeV. These values are compatible with those coming from the 
fit to the total mass spectrum. 

In addition to these waves, an extra wave, not required in the WA76 and WA91 analyses, 
is found to be necessary in the fit, namely the J p = 2~ a 2 (1320) 7r wave. The addition of this 
wave increases the log likelihood by 2700. The resulting J p = 2" a 2 (1320) 7r wave is shown in 
fig. |] where a broad enhancement is observed which peaks at 1.6 GeV. Other experiments have 
observed evidence for J PC = 2~ + states in this region. In the reaction 77 — > t/ti-k the Crystal 
Ball collaboration [17] reported a J PC = 2~ + state decaying to a2(1320)7r and ao(980)7r with 



a mass of 1881 ±32 ±40 MeV and a width of 221 ± 92 ± 44 MeV. In pp interactions the 
Crystal Barrel collaboration [18| reported the observation of two J PC = 2 h states in the r]7nr 



final state. The first state, called the 772(1620), had a mass of 1645 ± 14 ± 5 MeV and a width 
of 180 ± 40 ± 25 MeV and the second state, called the 772 (1875), had a mass of 1875 ± 20 ± 35 
MeV and a width of 250 ± 25 ± 45 MeV. 

As can be seen from fig. |2] the J p = 2~ a 2 (1320)7r wave observed in this experiment is 
consistent with the two r\ 2 resonances discussed above with both states decaying to a2(1320)7r. 
The masses and widths found from the fit are given in table [1|. 

In order to calculate the contribution of each resonance as a function of the dPr the distri- 
butions in fig. [5] have been fitted with the parameters of the resonances fixed to those obtained 
from the fits to the total data. The results of the fits are given in table 

As is observed from figs. [I] and |3] and from table ^ the /i(1285) signal almost disappears 
at small dPr- In addition, the 2 h a 2 (1320)7r signal is also suppressed at small dPx- This 
behaviour is consistent with the signals being due to standard qq states @. 

The ++ pp wave shows that the /o(1450) peak is still prominent at small dP?. If the peak 
is interpreted as the interference between the / (1300) and / (1500) then, as can be seen from 
table 0, both must have a similar dPx dependence. The broad peak at 2.0 GeV in the ++ pp 
wave, called the /o(2000), is less prominent at small dP? and is more consistent with other qq 
states. 

Both the J p = 2 + a 2 (1320)7r and J p = 2 + / 2 (1270)(7r7r),5 U , a , ;e distributions have a very 
similar dPj- behaviour consistent with there being just one resonance at 1.9 GeV with two 
decay modes. In fig. ^| the sum of these two waves is presented. As can be seen from this figure 
and from table ^| at small dPx the / 2 (1900) signal is still important. This is the first evidence 
of a non-zero spin resonance produced at small dPr and hence shows that the dPx effect is not 
just a J p = + filter. 

It is interesting to note that the prominence of the /o(1450) and /2(1900) signals in fig. |l]b) 
for dPx < 0.2 GeV is not only due to the fact that the / (1450) and / 2 (1900) signals are strong 
at small dPx but also due to the fact that the r/ 2 signals are suppressed. 

In summary, in the centrally produced 7r + 7r~7r + 7r~ mass spectrum two interesting structures 
are observed at 1.45 and 1.9 GeV. A spin analysis shows that the underlying resonance structure 
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is more complex. The peak at 1.45 GeV is found to have J PC = ++ and to decay to pp; it 
can either be described as a single Breit-Wigner or as an interference between the /o(1300) and 
/o(1500). The peak observed at 1.9 GeV appears as a broad enhancement in the mass spectrum 
and is found to decay to a 2 (1320)7r and / 2 (1270)7r7r with J PC = 2 ++ . Due to the difficulty in 
describing the individual a 2 (1320)7r and / 2 (1270)7T7r contributions by a Breit Wigner it is likely 
that this is one state, called the / 2 (1900), with two decay modes. Confirmation is found for 
two J = 2 _+ resonances, called the ?7 2 (1620) and 772 (1875) decaying to a 2 (1320)7r. There is 
also evidence for a broad state with J PC = ++ decaying to pp, called the / (2000). 
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Tables 



Table 1: Parameters of resonances in the fit to the n + n n + 7i mass spectrum and waves. 





Mass (MeV) 


Width (MeV) 


Observed 
decay mode 


I(J PC ) 


h (1285) 


1281 ± 1 


24 ± 3 




0(1++) 


/o(1300) 


1290 ± 15 


290 ± 30 


pnn 


0(0++) 


/o(1500) 


1510 ± 20 


120 ± 35 




0(0++) 


/o(2000) 


2020 ± 35 


410 ± 50 


pun 


0(0++) 


/ 2 (1900) 


1960 ± 30 


460 ± 40 


a 2 (1320)7r 

/ 2 (1270)7T7T 


0(2++) 


772(1620) 


1620 ± 20 


180 ± 25 


a 2 (1320)7r 


0(2-+) 


r? 2 (1875) 


1840 ± 25 


200 ± 40 


a 2 (1320)7r 


0(2-+) 
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Table 2: Resonance production as a function of dPx expressed as a percentage of its total 
contribution. 





dP T <0.2 GeV 


0.2<dP T <0.5 GeV 


dP T >0.5 GeV 


/i(1285) 


5.8 ± 0.4 


49.3 ± 1.3 


44.8 ± 0.8 


/o(1300) 


28.4 ± 1.0 


43.2 ± 1.0 


28.4 ± 2.0 


/o(1500) 


30.8 ± 1.0 


50.5 ± 1.0 


18.8 ± 2.0 


/o(2000) 


12.2 ± 2.0 


53.1 ± 2.0 


34.7 ± 2.5 


/ 2 (1900) 


31.1 ± 2.0 


49.2 ± 2.0 


19.7 ± 1.0 


?7 2 (1620) 


2.4 ± 1.0 


50.6 ± 3.4 


46.8 ± 3.3 


^ 2 (1875) 


0.4 ± 0.4 


39.3 ± 6.0 


60.3 ± 6.0 
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Figures 

Figure 1: The 7r + 7r~7r + 7r~ effective mass spectrum a) for the total data with fit using 3 Breit- 
Wigners b) for dP T < 0.2 GeV, c) for 0.2 < dP T < 0.5 GeV and d) for dP T > 0.5 GeV. 

Figure 2: Results of the spin parity analysis. The superimposed curves are the resonance 
contributions coming from the fits described in the text. 

Figure 3: Results of the spin parity analysis as a function of dPx- 
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